Asymmetric two-component Fermi gas with unequal masses 
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We analyze the zero temperature phase diagram for an asymmetric two-component Fermi gas as a 
function of mass anisotropy and population imbalance. We identify regions corresponding to normal, 
or uniform/non-uniform superfluid phases, and discuss topological quantum phase transitions in the 
Bardeen-Cooper-Schrieffer (BCS), unitarity and Bose-Einstein condensation (BEC) limits. Lastly, 
we derive the zero temperature low frequency and long wavelength collective excitation spectrum, 
and recover the Bogoliubov relation for weakly interacting dilute bosons in the BEC limit. 
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The evolution from Bardeen-Cooper-Schrieffer (BCS) 
to Bose-Einstein condensation (BEC) is a very impor- 
tant topic of current research for the condensed mat- 
ter, nuclear, atomic and molecular physics communi- 
ties. Recent advances in atomic physics have allowed 
for the study of superfluid properties in symmetric two- 
component fermion superfiuids (equal mass and equal 
population) as a function of scattering length, where the 
theoreticahypredicted crossover from BCS to BEC was 
observed 

Since the population of each component as well as the 
interaction strength between two components are exper- 
imentally tunable, these knobs enabled the study of the 
BCS to BEC evolution in asymmetric two-component 
fermion superfiuids (equal mass but unequal popula- 
tion) H0. In contrast with the crossover physics found 
in the symmetric case, these experiments have demon- 
strated the existence of phase transitions between nor- 
mal and superfluid phases, as well as phase separation 
between superfluid (paired) and normal (excess) fermions 
as a function of population imbalance. 

Arguably one of the next frontiers of exploration 
in cold Fermi gases is the study of asymmetric two- 
component fermion superfluidity (unequal masses and 
equal or unequal population) in two-species fermion mix- 
tures from the BCS to the BEC limit. Earlier works on 
two-species fermion mixtures were limited to the BCS 
regime However, very recently, the evolution 

from BCS to BEC was preliminarily addressed in ho- 
mogenous systems as a function of population imbalance 
and scattering length for 6 Li and 40 K mixture In ad- 
dition, the superfluid phase diagram of trapped systems 
at unitarity was also analyzed as a function of population 
imbalance and mass anisotropy (1 0| . 

In this manuscript, we study the BCS to BEC evolu- 
tion of asymmetric two-component fermion superfiuids as 
a function of population imbalance and mass anisotropy. 
Our results for a homogeneous system are as follows. 
We analyze the zero temperature phase diagram for an 
asymmetric two-component Fermi g function of 

mass anisotropy and population imbalance. We identify 
regions corresponding to normal, and uniform or non- 



uniform superfluid phases, and discuss topological quan- 
tum phase transitions in the BCS, unitarity and BEC 
limits. Lastly, we derive the zero temperature low fre- 
quency and long wavelength collective excitation spec- 
trum, and recover the Bogoliubov relation for weakly in- 
teracting dilute bosons in the BEC limit. 

To describe a dilute asymmetric two-component Fermi 
gas in three dimensions, we start from the pseudo-spin 
singlet Hamiltonian {% = 1) 



k,(T 



k,k',q 



where the pseudo-spin a labels the hyperfine states rep- 
resented by the creation operator a k a , and = 



Here, £ ki<7 = e k , CT - fi a , where e ki(T 



k / (2m a ) is the energy and /j, a is the chemical potential 
of the fermions. Notice that, we allow for the fermions to 
have different masses m„ and different populations con- 
trolled by independent chemical potentials fi a . The at- 
tractive fermion-fermion interaction V(k, k') can be writ- 
ten in a separable form as V(k, k') = — ^r^IV where 
g > 0, and T k = 1 for the s-wave contact interaction 
considered in this manuscript. 

The gaussian effective action for H is S'gauss = 
S + (/3/2)E g At(g)F- 1 (?)A(g) I where q = (q,v e ) with 
bosonic Matsubara frequency vg = 2£tt//3. Here, (3 = 
l/T, At(g) is the order parameter fluctuation field, and 
the matrix F~ 1 (q) is the inverse fluctuation propagator. 
The saddle point action is 



So = P 



|Aq| 



J2 {/?(&,+ - E K+ ) 

k 

Hn F (-E^)]+ln[n F (-E Kl )}},(2) 



where E k (T = | A^ 2 ) 1 / 2 + s CT £ k ,_ is the quasipar- 

ticle energy when sj = 1 or the negative of the quasi- 
hole energy when s± = —1, and £7 k ,± = C^M ^ ^k,i)/2. 
Here, A k = Aor k is the order parameter, nj?(^ kl0 -) is 
the Fermi distribution and £k,± = (£ k .f ± £ k ,j,)/2 = 
k 2 /(2m±) — fi±, where m± = 2m-\mi/ '(?7ij, ± m-\) and 
M± = (MT i Notice that m + is twice the reduced 



2 



mass of the f and J, fermions, and that the equal mass 
case corresponds to \rri—\ — * oo. The fluctuation term in 
the action leads to a correction to the thermodynamic 
potential, which can be written as f2 



gauss 



fin + Q 



fluct 



with fi„ - So/ (3 and O fluct = (l/P)J2 q \iidet[F~\q)/(3}. 

The saddle point condition SSo/SA^ = leads to an 
equation for the order parameter 



-E 



|r k | 2 x 



(3) 



where Ak,± = (<-fk,T ^ ■^k,|)/2 with X^ i<7 = 
tanh(/5£k i(J /2). As usual, we eliminate g in favor 
of the scattering length olf via the relation 1/g = 
-m+y/(47raF)+Ek |r k | 2 /(2e k , + ), where e k ,± = (e k)T ± 
£k,i)/2- The order parameter equation has to be solved 
self-consistently with number equations N a = —dQ,/dp a 
which have two contributions N a = N 0:(T + A^fl ucti(T - N _ a 
is the saddle point number equation given by 



N, 



0,er 



an 

dp a 
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and -/Vfi uc t j(T = —dflfiuct/d^a- is the fluctua- 
tion contribution to N given by -/Vfl uct a = 
-(l/^E.^IdetF-i^/a^j/detF- 1 ^). 

In order to analyze the phase diagram at T = 0, we 
solve the saddle point self-consistency (order parameter 
and number) equations as a function of population im- 
balance P = N-/N+ and mass anisotropy m r — m^/mi. 
Here, N± — (Nf ± -^VjJ/2 corresponding to the following 
momenta k F ± = (k F ^ ± k F j)/2. In addition, we check 
the stability of saddle point solutions for the uniform su- 
perfluid phase using two criteria. The first criterion re- 
quires that the curvature (9 2 f2 /9A§ of the saddle point 
thermodynamic potential Slo with respect to the saddle 
point parameter Aq to be positive, where 



dA 2 
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Here, y k ,± = (^ k , T ± y Kl )/2 with y Ka = 
sech 2 (/3_Ek.cr/2). When 9 2 O /i9Aq is negative, the uni- 
form saddle point solution does not correspond to a min- 
imum of Oo, and a non- uniform superfluid phase is fa- 
vored. This criterion is related to the positive dcfmitcncs 
criterion on the compressibility matrix k with elements 
n a ,cr' = — d 2 VLp/ (dp a dp a '). The second criterion requires 
that the superfluid density 



k 



(6) 



to be positive. When po(T) is negative, a spontaneously 
generated gradient of the phase of the order parameter 
appears, leading to a non-uniform superfluid phase. 




FIG. 1: Phase diagram of P = (7V T - iV|)/(iV T + N L ) versus 
m r = mf/raj on the BCS side when a) 1/(A;f,+<if) = —1 and 
b) l/(k F ,+a F ) = -0.25. 



Based on these two criteria, we construct the P versus 
to,, phase diagram for five sets of interaction strengths: 
l/(fcir-|_ op) = —1 and —0.25 on the BCS side shown in 
Fig. l/(fcp i+ a_F) = at unitarity shown in Fig. [21 
and l/(k F ^ + a F ) — 0.25 and 1 on the BEC side shown in 
Fig. |3 In these diagrams, the f (!) label always corre- 
sponds to lighter (heavier) mass such that lighter (heav- 
ier) fermions are in excess when P > (P < 0). Notice 
that this choice spans all possible population imbalances 
and mass ratios. In Figs. 11131 we indicate the regions of 
normal (N), and non-uniform (NU) or uniform (U) su- 
perfluid phases. The normal phase is characterized by 
a vanishing order parameter (Aq — 0), while the uni- 
form superfluid phase is characterized by po(0) > and 
d 2 ilo/dA 2 > 0. The non-uniform superfluid phase is 
characterized by po(0) < and/or d 2 £lo/dA 2 l < 0, and 
it should be of the LOFF-type having one wavevector 
modulation only near the BCS limit, although closer to 
unitarity, we expect the non-uniform phase to be sub- 
stantially different from the LOFF phases having spa- 
tial modulation that would encompass several wavevec- 
tors. However, from numerical calculations, the first 
criterion seems to be dominant for all parameter space 
and the non-uniform superfluid phase is characterized by 
d 2 ilo/dAQ < 0, which indicates possibly phase separa- 
tion. 

We also identify a topological quantum phase transi- 
tion, which is shown as dotted lines in Figs lll-il These 
phases are characterized by the number of zeros of S k ,| 
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FIG. 2: Phase diagram of P = (JV T - ATj.)/(iV t + iV;) versus 
m r = m-f/mj, at the unitarity when l/(fcj?+a.F) = 0. 

and -Ek.i (zero energy surfaces in momentum space) such 
that I) i?k,cr has no zeros and i?k,-o- has only one, and II) 
-Ek.a- has no zeros and -Ek,-<r has two zeros. Phase I (II) 
always appears to the left (right) of the dotted lines. No- 
tice that, the P = limit corresponds to case III, where 
i?k,o- has no zeros and is always positive. The transition 
from case II to case I occurring at the dotted lines is 
quantum in nature, and signatures of it could still be ob- 
served at finite temperatures through the measurement 
of momentum distributions However, phase II seems 
to lie always in the NU region and may not be accessible 
experimentally. 
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FIG. 3: Phase diagram of P = (JV T - iV;)/(JV T + N t ) versus 
m r = mi/mi on the BEC side when a) l/(kp t+ ap) = 0.25 
and b) l/(fcp,+aF) = 1. 

As shown in Fig. we find a small region of uni- 
form superfluidity on the BCS side for 1/ (k,F,+aF) = — 1 



only when the mass anisotropy is small and the lighter 
fcrmion are in excess {P > 0). Thus, mixtures consisting 
of 6 Li and 40 K [m r « 0.15) or 6 Li and 87 Sr (m r w 0.07) 
are good candidates for future experiments. In the rest 
of the phase diagram, we find a quantum phase transi- 
tion from the non-uniform superfluid to the normal phase 
beyond a critical population imbalance for both positive 
and negative P. The phase space of uniform superfluid- 
ity expands while that of the normal phase shrinks with 
increasing interaction strength as shown in Fig. ^p. 

This general trend continues into the unitarity limit 
\\/{kF,+a.F) = 0] as shown in Fig. |3 Since this limit is 
theoretically important as well as experimentally acces- 
sible, we suggest that Fermi mixtures corresponding to 
< m r < 0.45 have phase diagrams which are qual- 
itatively different from those of 0.45 < m r < 1 as a 
function of P. Thus, mixtures consisting of 6 Li and 2 H 
(to,. « 0.33), 6 Li and 25 Mg (m r « 0.24) or 40 K and 87 Sr 
[m r ~ 0.64) are also good candidates for future exper- 
iments. Notice that, our results for the case of equal 
masses (m r = 1) are in close agreement with recent 
MIT experiments Q in a trap. At unitarity, our non- 
uniform superfluid to normal state boundary occurs at 
P « ±0.73, and the MIT group obtains P » ±0.70(4) 
for their superfluid to normal boundary. 

Additional increase of interaction strength beyond uni- 
tarity on the BEC side leads to further expansion (shrink- 
age) of the uniform superfluid (normal) region as shown 
in Fig. |3| When heavier fermions are in excess (P < 
0), a uniform superfluid phase is not possible for any 
mass anisotropy until a critical interaction strength is 
reached. The critical interaction strength corresponds to 
1/ (kF,+o.F) ~ 0.8 for m r = 1. However, in the extreme 
BEC limit [l/(kp, + ap) » 1], only the uniform super- 
fluid phase exists even for P < (not shown). 

Next, we analyze gaussian fluctuations from which we 
extract the low frequency and long wavelength collective 
excitations at zero temperature. The collective excitation 
spectrum is determined by the poles of the propagator 
matrix F(g) determined by the condition det F^ 1 (q) = 0, 
when the usual analytic continuation ivi —> w + i0 + 
is performed. First, we express the fluctuation field as 
A(q) = [X(q) + iO(q)]/V2, where X(q) and 0(g) are am- 
plitude and phase fields, respectively. Then, we consider 
only the P = or kp,+ = kp^ = kp^ limit, and expand 
the matrix elements of F _1 (g) to second order in |q| and 
w to get 

, / A + C\q\ 2 -Dw 2 iBw \ 
^ w )=\ - lBw QH 2 -Rw 2 )■ W 

Thus, there are two branches for the collective excita- 
tions, but we focus only on the lowest energy one corre- 
ponding to the Goldstone mode with dispersion u>(q) = 
v |q|, where v = y/ AQ / {AR + B 2 ) is the speed of sound. 
Notice that extra care is required when P ^ since Lan- 
dau damping causes collective excitations to decay into 
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the two-quasiparticle continuum even for the s-wave case. 

The BCS limit is characterized by the criteria [i+ > 
and 6f,+ S> |Ao|. The expansion of the ma- 

trix elements to order |q| 2 and w 2 is performed un- 
der the condition [w, |q| 2 /(2m + )] <C |A |. The co- 
efhcient that couples phase and amplitude fields van- 
ish (B = 0) in this limit. Thus, there is no mix- 
ing between the phase and amplitude modes. The ze- 
roth order coefficient is A — T>. and the second order 
coefficients are C = Q/3 = Puf,tWf,i/(36| A | 2 ), and 
D = R/3 = £>/(12|A | 2 ). Here, v F<a = k F>a /m a is the 
Fermi velocity and T> = m + Vkp + / '(2tt 2 ) is the density 
of states per spin at the Fermi energy. Thus, we obtain 
v = \Z v F,1 v F,l/3 — vA^T' with v a = i'F,a/V3, which 
reduces to the Anderson-Bogoliubov relation when the 
masses are equal. 
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FIG. 4: Sound velocity v (in units of vf,+ ~ kF,+ /m+) versus 
m r = m^/mi for 1/(/cf,+cif) = —1 (solid line), 1/ (kF.+ap) = 
(solid squares) and l/(fc_F,+fflF) = 1 (hollow squares). Here, 
populations are equal (P — 0). 

On the other hand, the BEC limit is characterized by 
the criteria fi + < and £k.+ 3> |Ao|. The expansion of 
the matrix elements to order |q| 2 and w 2 is performed 
under the condition [w, |q| 2 /(2m + )] -C The coef- 

ficient B =^ indicates that the amplitude and phase 
fields are mixed. The zeroth order coefficient is A — 
k\ Ao| 2 /(2|/z+|), the first order coefficient is B = n, and 
the second order coefficients are C = Q = /c/[2(mj 
and D = R = rc/(8|/i + |), where k = V/(32^\fj, + \e F ,+)- 
Thus, we obtain v — |Ao|/a/4(to| + mjj|//+| = y/vfvl 
with v a = y / 27rn a aF /m 2 . Notice that the sound ve- 
locity is very small and its smallness is controlled by 
the scattering length ap. Furthermore, in the theory of 
weakly interacting dilute Bose gas, the sound velocity 
is given by vb — y 47ra bbtib/ ■ Making the iden- 
tification that the density of pairs is ns = the 
mass of the bound pairs is tub = m-f + mj, and that 
the Bose scattering length is cibb — {tub /m+)aF = 
[1 + TOj/(2mjJ + mi/(2mi)]aF, vb reduces to the well 
known Bogoliubov relation when the masses are equal. 



Therefore, the strongly interacting Fermi gas with two 
species can be described as a weakly interacting Bose 
gas at zero temperature as well as at finite tempera- 
tures 0- Notice that a B B reduces to a bb = 2a F for 
equal masses in the Born approximation, but a bet- 
ter estimate for as b can be found in the literature |l2( . 

In Fig. we show the sound velocity as a function of 
the mass ratio m r for three values of the scattering pa- 
rameter \/{kF,+a,F) = —1,0 and 1 corresponding to the 
BCS side [l/(kF,+a,F) = — 1], unitarity \l/(kF,+a F ) = 0], 
and to the BEC side [l/(k Fi +a F ) — 1]. Notice that 
the speed of sound could be measured for a given m r 
using similar techniques as in the single species case 
m r = 1 00. 

In summary, we analyzed the zero temperature phase 
diagram for an asymmetric two-component Fermi gas 
as a function of mass anisotropy and population imbal- 
ance. We identified regions corresponding to normal, and 
uniform or non-uniform superfluid phases, and discussed 
topological quantum phase transitions in the BCS, uni- 
tarity and BEC limits. Lastly, we derived the zero tem- 
perature low frequency and long wavelength collective 
excitation spectrum, and recovered the Bogoliubov re- 
lation for weakly interacting dilute bosons in the BEC 
limit. We thank NSF (DMR-0304380) for support. 

Note added: We would like to acknowledge valuable 
e-mail correspondence with S.-K. Yip from June to Au- 
gust of 2006 in connection with differences and similari- 
ties between the stability criteria for the phase diagrams 
of the original versions of references 0] and as well 
as of the present manuscript. In addition, recently, a few 
manuscripts have appeared addressing phase diagrams 
and their stability criteria for unequal population Fermi 
systems 0, |la, LLZl LLa, UJ) ■ These manuscripts have lead 
to a debate whether the compressibility criterion and the 
criterion based on the curvature of the thermodynamic 
potential with respect to the order parameter are equiv- 
alent or not. Our comments in connection to the current 
debate will be posted at a later date. 
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